Time-reversal symmetry breaking state near the surface of s-t-superconductor 
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The structure of superconducting order parameter near the surface of a two-band superconductor 
with s± order parameter in the bulk is theoretically investigated. The main parameter of the 
surface, which determines the appropriate physics is the coefficient of the interband scattering R12. 
For small R12 the superconducting order parameter is only suppressed to some extent near the 
surface for the both bands. For intermediate and strong interband scattering there are two possible 
non-trivial surface states of the order parameter: (i) purely real solution, where the symmetry of the 
superconducting state near the surface is changed from s± to conventional s++ and (ii) time-reversal 
symmetry breaking (TRB) state. In this state the order parameters in the two bands acquire phases 
</>i,2(a:) 7^ (0, 7r) upon approaching the surface. We argue that at low temperatures the TRB surface 
state can be more energetically favorable than the s± — > s++ time reversal symmetry conserving 
state (TR). For higher temperatures up to T c only the TR state can exist. The transition between 
the two temperature regions is rather sharp. Signatures of the transition between the TRB and 
the TR surface states can be detected by the measurements of the local density of states and the 
angle-resolved density of states. 

PACS numbers: 74.45.+C, 74.20.Rp, 74.70.Xa 



I. INTRODUCTION 



The discovery of a new family of iron-based high- 
temperature superconductors with distinct multi-orbital 
band structure^— has renewed interest to the problem 
of multi-band superconductivity, firstly discussed fifty 
years ago^£. In iron-pnictide^ and iron chalcogenides 7 - - — 
the Fermi surface (FS) has two hole-pockets around 
T = (0,0) and two electron-pockets around M — (tt, tt). 
The proximity of superconductivity to an antiferromag- 
netic phase suggests that the magnetic fluctuations play 
an important role in the understanding of pairing mecha- 
nism. It was propose d 11 ^ 2 that in this case the so-called 
s± superconducting state can be realized. It is charac- 
terized by the phase difference n between the supercon- 
ducting condensates arising on the hole Fermi surfaces 
around T point and the electron Fermi surfaces around 
M point. This state has been favored by a variety of 
models within random phase approximation (RPA)^~— 
and renormalization group techniques^—. 

The problem of inhomogeneities, such as disorder, sur- 
faces and interfaces, is substantially more complicated in 
these materials than in one-band systems and contains 
rich new physics, which is explored insufficiently by now. 
For simplicity we consider further only two Fermi surfaces 
(one electron FS and one hole FS). An incoming quasi- 
particle belonging to one Fermi surface (we also use term 
"band" below) can be scattered by the impurity (surface, 
interface) as into the same band (intraband scattering), 
so as into the other band (interband scattering). Inter- 
play of the interband scattering with the different signs of 
the superconducting order parameter (OP) on different 
Fermi surfaces can lead to various new physical phenom- 
ena. For example, it has been shown that in the s± state 
any fully interband nonmagnetic impurity suppresses T c 
in the same way as a magnetic impurity in a single band 



superconductor™. However, a much slower pairbreaking 
rate can be achieved if one assumes that the scattering is 
primarily intraband rather than interband. In the pure 
intraband scattering limit Anderson's theorem can be ap- 
plied and, consequently, no T c suppression occurs. The 
rate of T c suppression is, therefore, very sensitive to the 
relative strength of the intra- and interband scattering 
rates. Further, it was shown' 19 that depending on the 
ratio of inter- to intraband pairing constants there are 
two possible types of s± superconductivity. For the first 
one T c is suppressed as disorder is increased and vanishes 
at a critical value of the scattering rate. For the second 
type T c tends to a finite value upon increasing of dis- 
order, while the order parameter undergoes a transition 
from s± to a gapless state and then to a fully gapped 
,s ++ state. 

Surface and interface phenomena in s±- 
superconductors have also attracted considerable 
recent attention. The formation of bound states at a free 
surface of an s±-superconductor— ~— , at an S±/N22r— , 
an N/S/S± junction 28 and at Josephson junctions 
including 5±-superconductor o 27 ' 29 was investigated the- 
oretically. In particular, the finite energy subgap bound 
states were found and their influence on the conductance 
spectra and Josephson current was investigated. It 
was found that the particular energies of the bound 
states are very sensitive to the relative strengths of the 
intraband and interband scatterings. 

Further, it was found 30 that the spatial behavior of the 
OP at a surface of a s± superconductor can not be re- 
duced to a trivial suppression. If the interband scattering 
at a surface is of the order of the intraband one or domi- 
nates it, the symmetry of the superconducting state near 
the surface can be changed from s± to conventional 
The spatial region of existing of this surface conventional 
superconductivity is very sensitive to the relative values 
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of interband and intraband pairing potentials. The tran- 
sition between the bulk s± an surface s ++ superconduc- 
tivity regions occurs through a region of suppressed order 
parameter, to some extent analogously to the discussed 
above gapless superconductivity in the disorder problem. 

In the present paper we continue investigation of the 
OP spatial behavior near a surface of a s± supercon- 
ductor. It is found that besides the discribed above 
purely real solution for the order parameter there is an- 
other solution, which starts from the same s± state in 
the bulk, but manifests qualitatevely different behavior 
near the surface. This state is time-reversal symmetry 
breaking (TRB): upon approaching the surface the or- 
der parameters in the two bands (called 1 and 2) acquire 
phases 4>\.2(x) ^ (0, it), while the bulk values <pi = and 
<t>2 = 7r restore within a few superconducting coherence 
lengths from the surface. We argue that this TRB sur- 
face state can be more energetically favorable than the 
discussed above s± — > s ++ time reversal symmetry con- 
serving state (TR). The TRB surface state only occurs 
at low temperatures, while for higher temperatures up to 
T c only the TR state can exist. The transition between 
the two temperature regions is rather sharp. Further, we 
show that the local density of states (LDOS) spectra are 
very sensitive to the particular surface state of the or- 
der parameter and, therefore, experience a sharp change 
when temperature crosses the transition point. There- 
fore, the transition between TRB and TR surface states 
can be detected by investigation the LDOS. 

The possibility of TRB state has been already dis- 
cussed in literature 3 - 1 - - — . However, the physical mecha- 
nism, underlying this TRB state is a frustration between 
three or more Josephson-coupled bands. On the contrary, 
we consider only two Fermi pockets and the surface TRB 
state is generated by the interband scattering between 
them. 

The paper is organized as follows. The model sys- 
tem we study and our theoretical method are described 
in Sec. HU In Sec. IIIII we represent the results of the 
self-consistent calculations of the superconducting order 
parameter for our surface problem and discuss the TRB 
state. In Sec. IIVI the LDOS is calculated and its depen- 
dence on temperature is analyzed. Our conclusions are 
given in Sec. [Vj 



II. MODEL AND METHOD 

We consider an impenetrable surface of a clean two- 
band superconductor. The OP is assumed to be of s±- 
symmetry in the bulk of the superconductor, that is the 
phase difference between the OP's in the two bands is 
7r. It is supposed that an incoming quasiparticle from 
band 1,2 can be scattered by the surface as into the same 
band (intraband scattering), so as into the other band 
(interband scattering). In order to maximize the effects 
generated by the surface scattering we assume that the 
c-axis is parallel to the surface. 



We make use of the quasiclassical theory of supercon- 
ductivity, where all the relevant physical information is 
contained in the quasiclassical Green function f)i(e,pf, x) 
for a given quasiparticle trajectory. Here e is the quasi- 
particle energy measured from the chemical potential, pf 
is the momentum on the Fermi surface (that can have 
several branches), corresponding to the considered tra- 
jectory, x is the spatial coordinate along the normal to 
the surface and i = 1, 2 is the band index. Quasiclassical 
Green function is a 2 x 2 matrix in particle-hole space, 
that is denoted by the symbol The equation of mo- 
tion for gi(e,pf, x) is the Eilenberger equation subject to 
the normalization conditio n 43 ' 44 . For superconductivity 
of s±-type electrons from different bands cannot form a 
pair. Therefore, the Eilenberger equations corresponding 
to the bands 1 and 2 are independent. The trajectories 
belonging to the different bands can only be entangled 
by the surface, which enters the quasiclassical theory in 
the form of effective boundary conditions connecting the 
incident and outgoing trajectories. 

However, owing to the normalization condition for the 
quasiclassical propagator, the boundary conditions for 
the quasiclassical Green functions are formulated as non- 
linear equations^ 5 - - — . For the multiband case, consid- 
ered here, they are rather complicated and their practi- 
cal use is limited. For this reason in the present work 
we make use of the quasiclassical formalism in terms 
of so-called Riccati amplitude a 48 ' 49 , which allows an ex- 
plicit formulation of boundary conditions^ 9 - - — . The re- 
tarded Green function gi(e,pf,x), which is enough for 
a complete description of an equilibrium system, can be 
parametrized via two Riccati amplitudes (coherence func- 
tions) ji(e,pf,x) and 7,;(e, Pf , x) (in the present paper 
we follow the notations of Refs. |49| and (53|). The co- 
herence functions obey the Riccati-type transport equa- 
tions. In the considered here case of two-band clean s±- 
superconductor the equations for the two bands are in- 
dependent and read as follows 

iv ix d x ji + 2e-fi = -A*(^)7 4 2 - , (1) 

li(e,Pf,x) = i*{-e,-pf,x) . (2) 

Here Vi x is the normal to the surface Fermi velocity com- 
ponent for the quasiparticle belonging to band i. A.i(x) 
stands for the OP in the i-th band and should be found 
self-consistently. 

Let us suppose that the surface is located at x = 
and the superconductor occupies the halfspace x > 0. 
For the sake of simplicity we assume that the surface 
is atomically clean and, consequently, conserves paral- 
lel momentum component. Then there are four quasi- 
particle trajectories, which are involved in each surface 
scattering event. These are two incoming trajectories be- 
longing to the bands 1,2 (with Vi X < 0) and two outgoing 
ones (with Vi X > 0). It can be show n 49 ' 53 that the coher- 
ence function ji(e,pf,x), corresponding to the incoming 
trajectory can be unambiguously calculated making use 
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of Eq. (fTJ) and starting from its asymptotic value in the 
bulk 



Tip = -RoTifc + -Ri2T2fc- 



7» 



A^sgne 



\e\ + y/(e + i5)*-Af ' 



(3) 



where A\ is the bulk value of the OP in the appropriate 
band, 5 > is an infinitesimal. The coherence function 
7i(e,p/,x) is determined unambiguously by the asymp- 
totic conditions for the outgoing trajectories and can be 
obtained according to Eqs. ([I]),©. 

Otherwise, the coherence functions ji(s,pf,x) for the 
outgoing trajectories and, correspondingly, ji(e,pf,x) 
for the incoming ones should be calculated from Eq. (fTJ) 
supplemented by the boundary conditions at the surface 
and Eq. ([2j . The surface is described by the normal state 
scattering matrix for particle- like excitations, denoted by 
S and for hole-like excitations, denoted by S. The scat- 
tering matrix S have elements Sk iPj , which connect out- 
going quasiparticles from band i with momentum ki to 
the incoming ones belonging to band j with momentum 
Pj. Here and below all the momenta corresponding to 
the incoming trajectories are denoted by letter p and all 
the momenta for the outgoing quasiparticles are denoted 
by k. For the model we consider S is a 2 x 2-matrix in the 
trajectory space (for the particular value of the momen- 
tum component parallel to the surface, which is denoted 
by It obeys the unitary condition SS^ = 1. We also 
assume that the surface hamiltonian by itself posesses 
the time-reversal symmetry, which imposes an additional 
constraint on the scattering matrix S(pn) = S tr (— P||) 47 - 
In the absence of spin-orbit interaction the i'-matrix el- 
ements are only functions of |pn|- Then, without loss of 
generality the scattering matrix can be parameterized by 
three quantities R12, O and a as follows 
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where Rq and R12 are coefficients of intraband and inter- 
band reflection, respectively They obey the constraint 
Rq + R12 = 1. The sign factor a = ±1 and the phase 
factor O do not enter the boundary conditions and, conse- 
quently, do not influence all the obtained results. While 
in general the scattering matrix elements are functions 
of |p|||, we disregard this dependence in order to simplify 

the analysis. The scattering matrix S for hole-like excita- 
tions is connected to S by the relation S(pu ) = S tr (— pn), 

that is in the case we consider S — S. 

From the general boundary conditions^, which are 
also valid for a multiband system, one can obtain the 
explicit values of the coherence functions 7j(e,fe,a; = 0) 
and Ti( e : Pi x — 0) via the scattering matrix elements and 
the values of the coherence functions ji(e,p,x — 0) and 
7i(e, fe, x = 0) at the surface. They read as follows 




Tife = ^07ip+-Ri2T2p- 



-Rp^l272fc(7lp - T2 P ) 2 
1 + T2fc (-Rl2Tlp + -Ro72p) 



flo-Rl272p(7lfc - 72fc) 2 
1 + 72p (i?127lfc + R<)l2k) 



(6) 



Here the arguments (e, x = 0) of all the coherence func- 
tions are omitted for brevity, 7^ = 7* (p) and 7, p = 7, (p) 
and the analogous notations are used for ji(k) and Ti(fe)- 
Quantities jip and 7ik, entering Eqs. (JSJ) and (J6j> are to be 
calculated from Eqs. (flj and @ supplemented by the ap- 
propriate asymptotic condition. The coherence functions 
72fc and 72 P are obtained by the interchanging 1 -f-> 2 in 
all the coherence function band indices at the right-hand 
side of Eqs. ([5]) and ©, respectively. 

Now, substituting the coherence functions into the self- 
consistency equation 



1 + Tip/ Tip/ / 



(7) 



p/ 



, (5) 



we iterate system (QJ-([3]), ©-(JT]) until it converges. In 
Eq. (J7J Ajj < is the dimensionless pairing potential for 
band i and A12, A 2 i are the dimensionless interband pair- 
scattering potential. For simplicity we assume that the 
normal state densities of state coincide for the two FS: 
Nf,i — Nf,2- In this case A12 = A21. We choose A12 > 0, 
which stabilizes s± OP in the bulk. The Matsubara fre- 
quencies e n enter the coherence functions via the substi- 
tution e + iS — > ie n . (...) P j: means the anomalous Green 
function averaged over the entire Fermi surface, that is 
Pf incorporates as the incoming trajectories p, so as the 
outgoing ones k. For concreteness we suppose the Fermi 
surface to be cylindrical for the each band. However, 
our results do not qualitatively sensitive to the described 
above simplifying assumptions. 



III. SELF-CONSISTENT ORDER PARAMETER: 
TRB STATE 

The spatial profiles of the OP calculated according to 
the described above technique at T — 0.34T C , are rep- 
resented in Fig. [T] We assume that in the bulk \A\\ > 
|A|| and look for the solutions in the form Ai^a;) = 
|Ai ) 2(a;)| exp[i0i j2 (a;)]. Panel (a) of Fig. [1] demonstrates 
the spatial profiles of OP absolute values for bands 1 and 
2, while the corresponding phases are plotted in panel 
(b). It is seen from the Figure that at this temperature 
the self-consistency equation has two different solutions. 
One of them is the purely real (TR) solution, where the 
transition from s± to s ++ state at a distance ~ 0.25£i 
takes place. Here £1 = v\J A\(T = 0) is the supercon- 
ducting coherence length for band 1. This solution is 
characterized by going to zero ^2(2;) | and by abrupt 
jump of the phase 4>2{x) from 7r to 0. The other solution 
is time- reversal symmetry breaking. |Aa(a;)| does not ap- 
proach zero, while Ai^ acquire phases </>i,2(a;) 7^ (0,tt) 
in some region near the surface. The conjugated solu- 
tion characterized by — 4>x,2(x) is always exists, but is 
not represented in Fig. [1] 
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Due to the presence of the spatially dependent OP 
phase the TRB state is accompanied by the currents, 
flowing perpendicular to the surface and transferred by 
electrons from the both bands. However, the currents, 
transferred by the electrons from band 1 and from band 
2 exactly compensate each other at any distance x from 
the surface. As a result, the total current, generated by 
the spatially dependent OP phase in the TRB state, is 
zero, as it should be due to the charge conservation. 




0.5 1 1.5 2 2.5 

FIG. 1. The OP profiles as functions of the spatial coordi- 
nate x at T = 0.34. Panel (a) shows the OP absolute val- 
ues: \Aj RB \ (black solid line); |A;P S | (dashed); \Aj R \ (gray 
solid); |Aa fl | (dotted). The OP phases are demonstrated in 
panel (b): </>f fls (black solid); ^ RB (dashed); (gray 
solid); (j>2 R (dotted). R 12 = 0.5. In this Fig. and be- 
low throughout the paper An = 0.2771, A22 = 0.2241, and 
A12 = —0.004; all the spatial coordinates are measured in 
units of £1; all the energies, including the OP absolute val- 
ues, are measured in units of A\(T = 0) and temperature is 
measured in units of T c . 

The temperature evolution of the TRB solution is rep- 
resented in Figs. [2] and [3] Fig. [5] shows the OP absolute 
values, and Fig. [3] demonstrates the evolution of the OP 
phases. It is seen from the Figures that there is a rather 
narrow temperature region from w 0.3T C to w 0.4T C , 
where the OP profiles change qualitatively from the TRB 
state with nonzero phase to the real TR solution. The 
TRB state completely disappears above T ~ 0.42T C . At 
the same time the TR solution does not manifest any 
qualitative changes in the whole temperature region from 
to T c . 

In order to more explicitly show the temperature tran- 
sition between the TRB and the TR states in Fig. 2] we 
plot the OP phases <pi t 2 for the TRB state at the sur- 
face x = as functions of temperature. The dashed and 
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FIG. 2. The OP absolute values for the TRB state as func- 
tions of the spatial coordinate x. Panel (a) shows |Ai| and 
panel (b) represents |A2|. For the both panels different curves 
correspond to different temperatures: T = 0.30 (black solid), 
T = 0.34 (dashed), T = 0.38 (dashed-dotted) and T = 0.42 
(gray solid). R12 = 0.5. 
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FIG. 3. The OP phases <j>i (a) and <j>2 (b) as functions of 
the spatial coordinate x. The temperatures and the other 
parameters are the same as in Fig. [2] 



dotted lines represent phases 4>\{x — 0) and = 0), 
respectively, calculated for R12 = 0.5. The black solid 
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line corresponds to the limit of the full interband scat- 
tering R 12 = 1. In this case the values of 4>\{x = 0) 
and <p2 (x = 0) approximately coincide. The sharp con- 
version of the TRB surface state to the TR s ++ sur- 
face state is clearly seen at T » 0.4T C . It is important 
that we consider the case of weak interband coupling 
| A12 ] "C | An j j | A22 1 ■ Such a choice of parameters is con- 
sistent with the experimental estimates of the coupling 
constants for FeSe^. For more strong interband coupling 
the temperature transition between the TRB and the TR 
states would be more smeared. 




0.2 0.4 0.6 0.8 1 

T/T c 

FIG. 4. OP phases <j!>i,2 for the TRB state at the surface 
x — as functions of temperature. Dashed line: <fi\(x = 0) 
for R12 = 0.5, dotted line: <\>i(x — 0) for R12 = 0.5 and black 
solid line: 4>i(x = Q) = 4>2(% = 0) for R12 = 1. 

For the regime of weak interband coupling there are 
simple physical arguments indicating that the TRB state 
at low temperatures is more energetically favorable than 
the TR one. Let us consider the spatial region, where the 
OP A 2 (x) rea ', corresponding to the TR solution crosses 
zero point. It is the region giving the main difference 
between the free energies of the TRB and the TR states. 
Let us imagine that A2(x) real acquires small constant 
imaginary part iSA2 in the vicinity of the zero crossing 
point. Then, it is obvious that the leading (quadratic) 
term of GL free energy F ~ a%\ Ai| 2 +a 2 |A 2 | 2 + &(A^A 2 + 
AiA|) decreases for T < T c2 due to the negative contri- 
bution of the term a 2 |A 2 | 2 and increases for T > T c2 due 
to the positive contribution of this term, while the cou- 
pling term remains unaffected by the small imaginary 
part. Here T C 2 is the bulk critical temperature of the 
lesser order parameter A 2 under the condition of no inter- 
band coupling A 12 = (See Fig. [5]). Therefore, this rough 
consideration indicates that the TR state is unstable at 
low temperatures. It is the TRB state that becomes en- 
ergetically favorable at T < T c2 . Indeed, from Fig. [4] it 
is seen that the TRB state disappears at T « T c2 . For 
the regime of strong interband coupling all the above ar- 
guments are hardly applicable and in order to conclude 
which state (TRB or TR) is more favorable the more 
rigorous consideration of free energy functional is neces- 
sary. However, we do not consider the strong interband 
coupling regime in the paper. 




0.2 0.4 0.6 0.8 1 

T/T c 

FIG. 5. Bulk values Ai and A§ as functions of temperature 
for the system under consideration. The dashed lines show Ai 
and A2 under the condition of no interband coupling A12 = 0. 



IV. LDOS 

Now we discuss how the described above OP spatial be- 
havior affects the LDOS near the surface. LDOS p(x, e) 
is calculated via the coherence functions as follows 

i hp/ hp/ I pf 

As it was already discussed in the literature, if i?i 2 7^ 
there are surface bound states in the system, which man- 
ifest themselves as well-pronounced peaks in the LDOS. 
We begin by discussing the case of pure interband scatter- 
ing i?i 2 — > 1. Although this limit can hardly be realized 
in real experiments, the influence of the surface OP state 
on the LDOS is the most striking here. Then we turn 
to discussion of the effects taking place at more realistic 
intermediate values of the interband scattering. 

It was found^ that for the TR solution at i?i 2 — s- 1 
the bound state energy tends to zero. For this case the 
LDOS is dominated by very strong zero-energy peak. If 
the TRB surface state is realized in the system, then this 
is not the case. Instead of well-pronounced zero-energy 
peak the LDOS is characterized by two symmetric finitc- 
energy bound state peaks, which are very close to the 
smaller gap edges. This remarkable fact can be used in 
experiment in order to answer the question if there is the 
temperature transition between the TRB and the TR 
s ++ surface states in the system. Indeed, if the TRB is 
realized at low temperature, then at some intermediate 
temperature it should evolve rather sharply to the TR 
s+4- surface state. This transition is accompanied by the 
corresponding sharp modification of the LDOS spectra: 
two symmetric finite energy peaks divided by U-shaped 
inner gap evolve into the zero-bias peak. This is demon- 
trated in Fig. [6] It is worth noting that the main evolu- 
tion of the LDOS occurs in a rather narrow temperature 
interval, where the transition between the TRB and the 
,s ++ surface states takes place. The transition region is 
represented in panels (b)-(e) of Fig. [6l Beyond this tem- 
perature interval [panels (a) and (f)] the modification of 
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FIG. 6. LDOS at the surface x = as a function of quasiparticle energy e, which is measured in units of Aj(T = 0). The case 
of full interband scattering R12 = 1. 



the LDOS is not qualitative and mainly connected to the 
temperature evolution of the bulk superconducting gaps. 
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FIG. 7. Angle-resolved density of states at x — for the 
TRB state as a function of quasiparticle energy e. For all the 
panels: different curves correspond to different values of the 
angle 6 between the quasiparticle trajectory and the normal 
to the surface. 9 increases from bottom to top. (a) incoming 
trajectories for band 1; (b) outgoing trajectories for band 1; 
(c) incoming trajectories for band 2; (d) outgoing trajectories 
for band 2. Full interband scattering R12 = 1, T = 0.34. 

The LDOS behavior typical for the TRB state and 
demonstrated in panels (a)-(c) of Fig. [6] is quite simi- 
lar to the LDOS spectra corresponding to the ordinary 
s± surface state, where the bulk s± OP is only sup- 
pressed to some extent near the surface. This surface 
state should be realized for weak enough interband scat- 
tering i?i2 < 0.2 30 . Having in mind that it is very hard to 
estimate the value of the interband scattering R12 for the 
particular experimental situation the question emerges: 
is it possible to conclude on the basis of the LDOS spectra 
which surface state (suppressed s± or TRB) takes place 
in the system? The answer can be obtained from obser- 
vation the temperature evolution of the LDOS: while for 



the TRB state the discussed above sharp modification of 
the LDOS should be observed, there should be no quali- 
tative modification of the LDOS spectra upon increasing 
of temperature for s± surface state. 

The second signature of the TRB state can be found 
in the angle-resolved density of states and is closely con- 
nected to the time-reversal symmetry breaking nature of 
this state. In Fig. [7] we demonstrate the angle-resolved 
density of states for the case of full interband scattering. 
Panels (a) and (b) correspond to incoming and outgo- 
ing trajectories, respectively, in band 1, and panels (c) 
and (d) represent incoming and outgoing trajectories for 
band 2. It is clearly seen that the spectra are asymmet- 
ric with respect to zero energy. The reason is the follow- 
ing. The full interband scattering is a special case when 
not all the four trajectories (incoming in bands 1 and 2 
and outgoing in bands 1 and 2) are involved in the same 
scattering event at the surface, but the incoming trajec- 
tory from band 1 is only entangled with the outgoing one 
from band 2 (the first closed loop) and vice versa: the 
incoming trajectory from band 2 is only entangled with 
the outgoing one from band 1 (the second closed loop). 
These two closed loops are connected by time-reversal 
transformation, if one disregards the difference between 
p\\ and —p\\ (there are no reasons for such a difference 
because we do not study any effects of spin-orbit interac- 
tion and so on). The energies of the bound states, living 
on these loops, are connected by sign change. This is a 
manifestation of the TRB nature of the surface state. It 
is worth noting here that the substitution ^2 — > —01,2) 
that is the choice of the conjugated OP solution, leads 
to the sign change of all the bound state peaks in Fig. [7] 
However, integration over the Fermi surface results in 
disappearance of the difference between the two possible 
OP solutions in the LDOS. 

Now we turn to the more realistic case of the inter- 
mediate interband scattering R12 = 0.5. For the inter- 
mediate value of interband scattering there are finite en- 
ergy bound states in the system, which merge with the 



7 





j ■ 



- ' 




A. 










































1 1 




























FIG. 9. The same as in Fig. [Jj but for the intermediate inter- 
band scattering R12 = 0.5. 



smaller gap edges when R 12 — > 0. The LDOS behavior for 
the low-temperature TRB state and its evolution upon 
transition to the TR s ++ surface state is demonstrated 
in Fig. [8] Although the zero-energy peak is absent for 
R12 < 1, the qualitative difference between the LDOS 
spectra in the TRB and the s++ state exists. While in 
the TRB state the bound state peaks are divided by the 
clearly defined U-shaped gap, the transition to the s ++ 
state results in transforming this U-shaped gap into V- 
shaped behavior, which is known to be more typical for 
the superconductors with OP nodes at the Fermi surface. 
As it was shown—i, this fact is a consequence of spatial 
line of OP nodes appearing if the OP sign reversal takes 
place near the surface. Exactly as it was for the case of 
full interband scattering, the main qualitative evolution 
of the LDOS occurs in a rather narrow temperature in- 
terval, where the transition between the TRB and the 
s ++ surface states takes place. 

In this case angle-resolved spectra, calculated for the 
TRB state, also manifest the energy asymmetry. How- 
ever, in contrast to the limit of full interband scattering, 



now the both bound states (corresponding to positive 
and negative energies) are present for all the trajectories, 
but the corresponding peaks in the density of states have 
asymmetric heights. This is a consequence of the fact 
that for the intermediate interband scattering all four 
trajectories are involved in each scattering event at the 
surface. 

It is seen from Figs. [7J and ® that the bound states 
are practically dispersionless for R\2 = 1 and are only 
slightly dispersive for R12 = 0.5. This is a consequence of 
our simplified model for the surface: we assume that the 
scattering matrix elements do not depend on the direc- 
tion of of the incoming quasiparticle momentum. In order 
to check the generality of our main conclusions we have 
considered a number of momentum-dependent models for 
the scattering matrix. We have found that, although this 
leads to more dispersive angle-resolved spectra, the main 
qualitative effects, concerning the sharp modification of 
the LDOS upon varying temperature and the asymmetric 
angle-resolved spectra in the TRB state, are not sensitive 
to the details of the particular model. 



V. CONCLUSIONS 

On the basis of the self-consistent quasiclassical ap- 
proach we have theoretically investigated which types of 
superconducting order parameters are possible near the 
surface of a two-band superconductor with s± order pa- 
rameter in the bulk. It is found that the answer is mainly 
determined by the value of interband surface scattering 
with respest to the intraband one. If the intraband sur- 
face scattering Rq dominates the interband one R12, that 
is Rq is close to unity, then the superconducting OP is 
only suppressed to some extent near the surface for the 
both bands. For intermediate and strong interband scat- 
tering there are two possible non-trivial surface states of 
the order parameter. Further the particular answer de- 
pends on the ratio between the intraband and interband 
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pairing constants. We have focused on the case of weak 
interband coupling. It is found that at low temperatures 
the self-consistency equation has two solutions: (i) purely 
real solution, where the symmetry of the superconducting 
state near the surface is changed from s± to conventional 
s++. In this state the OP belonging to band 2 (with the 
smaller absolute value of the bulk OP) changes its sign 
near the surface via crossing zero value at some distance 
from the surface, (ii) Time-reversal symmetry breaking 
solution: upon approaching the surface the order param- 
eters in the two bands acquire phases fa^fa) 7^ (0,7r), 
while the bulk values fa — and fa = tt restore within 
a few superconducting coherence lengths from the sur- 
face. We argue that at low temperatures the TRB sur- 
face state can be more energetically favorable than the 
s± — > ,s ++ time reversal symmetry conserving state. For 
higher temperatures up to T c only the TR state can ex- 
ist. The transition between the two temperature regions 
is rather sharp. 

Further, we show that the local density of states is very 
sensitive to the particular surface state of the order pa- 
rameter and, therefore, experience a sharp change upon 



varying temperature. For the case of full interband scat- 
tering the finite-energy peaks, which are located close to 
the smaller gap edges and are typical for the TRB state, 
are changed by the zero-energy peak upon increasing of 
the temperature. For the case of intermediate interband 
scattering the zero-energy peak cannot exist, but the U- 
shaped low-energy LDOS behavior is changed by the V- 
shaped LDOS behavior at the transition from the TRB 
to the s++ surface state. In addition, the angle-resolved 
spectra manifest the energy asymmetry in the TRB state. 
Therefore, the signatures of the transition between the 
TRB and the s++ surface states can be detected by the 
density of states measurements. 
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